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ABSTRACT 

We investigate the potential of the Square Kilometre Array (SKA) for measuring the mag- 
netic fields in clusters of galaxies via Faraday rotation of background polarised sources. The 
populations of clusters and radio sources are derived from an analytical cosmological model, 
combined with an extrapolation of current observational constraints. We adopt an empirical 
model for the Faraday screen in individual clusters, gauged to observations of nearby clusters 
and extrapolate the polarisation properties for the radio source population from the NRAO 
VLA Sky Survey. We find that about 10 per cent of the sky is covered by a significant ex- 
tragalactic Faraday screen. Most of it has rotation measures between 10 and 100 rad mT^. 
We argue that the cluster centres should have up to about 5000 rad vcT^. We show that the 
proposed mid frequency aperture array of the SKA as well as the lowest band of the SKA 
dish array are well suited to make measurements for most of these rotation measure values, 
typically requiring a signal-to-noise of ten. We calculate the spacing of sources forming a grid 
for the purpose of measuring foreground rotation measures: it reaches a spacing of 36 arcsec 
for a 100 hour SKA observation per field. We also calculate the statistics for background RM 
measurements in clusters of galaxies. We find that a first phase of the SKA would allow us 
to take stacking experiments out to high redshifts (> 1), and provide improved magnetic field 
structure measurements for individual nearby clusters. The full SKA aperture array would be 
able to make very detailed magnetic field structure measurements of clusters with more than 
100 background sources per cluster up to a redshift of 0.5 and more than 1000 background 
sources per cluster for nearby clusters, and could for reasonable assumptions about future 
measurements of electron densities in high redshift clusters constrain the power law index for 
the magnetic field evolution to better than Am — 0.4, if the magnetic field in clusters should 

follow Boc 
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1 INTRODUCTION 

The plasma between the galaxies in clusters is magnetised. Ob- 
servations at radio frequencies are our main quantitative probe 
of cluster magnetic fields. This is done by three main meth- 
ods: equipartition argument s applied to cluster halo radio sources 
(e.g. iGiovannini et al.lll993h . comparison of X-ray Inverse Comp- 
ton (IC) with the radio synchrotron data in radio sources (e.g. 
iBagchi et aljl99^ . and Faraday rotation of polarised radio sources 
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within or behind the cluster, combined with X-ray data to deter - 



mine the electron densities jKim et alj[T99ll; Iciarke et al 



The results have been reviewed reg ularlv ( Carilli & Tavloij |2002l : 



2001) 



Govoni & Feretti 2004; Feretti & Giovannini 2008). Galaxy clus- 
ters are found to possess magnetic fields of typically a few, up to of 
order 10 /iG. 



The structure within cluster magnetic fields has been de- 
duced from Faraday r otation of extended c luster radio sources 
bet & EnfilinI l2005l : iGovoni et alj 12004 
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|2008l). For the few cases observed so far, the data is 



consistent with a Kolmogorov-type power spectrum. Because the 
field is frozen into the plasma in the high conductivity limit, the 
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Table 1. Comparison of radio telescopes — derived properties. 



Name 


VLA" 


SKA AA SKA Dishes 


SKA Dishes Phase 1 


Acff/Tiys'' 


70-180 


10,000 


10,000 


1,200 




0.3-0.34, 1.24-1.7, 








frequencies of interest / GHz 


4.5-5, 8.1-8.8 


0.3-1 


0.8-10 


0.8-10 


instantaneous bandwidth / MHz 


86 


700 


250 


250 


l-li sensitivity / fxiy 


17 


0.18 


0.29 


2.5 


field of view / square degree*^ 


{0.7//)2 


250 


2 


20'' 



^ See:lhttp : //mm . via . nrao . edu/astro/guld es/vlas/current/ 1 

Effective collecting area over system temperature. 

For the dish telescopes, given at the bottom of the frequency range; / denotes the observing frequency in GHz. 
Assuming phased an'ay feeds. 



Structure of the field directly traces the underlying kinematics of 
the gas. 

Current studies of magnetic fields in clusters are limited to a 
few nearby objects (z < 0.1). Hence, we have no means to answer 
questions regarding the origin and cosmological evolution of these 
fields. Correspondingly, few cosmological structure formation sim- 
ulations have included magnetic fields. To date, none have the res- 
olution to address directly the field evolution in galaxy clusters. 
However, from the eddy turnover times, and the high efficiency of 
the turbulent dynamo, one might expect that the magnetic energy 
should always be in equipartition with the turbulent kinetic energy 
(Ryu, private communication). Hence, studying the cosmological 
evolution of fields in clusters might not only lead to the origin of 
magnetism, but also trace cluster kinematics. 

Magnetic fields in clusters of galaxies might be connected 
to the cooling flow problem. If cooling flows in nearby objects 
are quenched by weak jets, as has been often suggested (e.g. 
iKaiser & Binnevll2003h . the magnetic fields will be amplified up 
to equipartition with the enhanced kinetic energy level provided by 
the very same jets. 

The Faraday rotation method is expected to benefit substan- 
tially from new developments of radio telescopes over the com- 
ing decade. The Australian ASKAP project, as well as the South 
African MEERKAT project are expected to start operation within 
a few years both instruments will have sensitivities competitive 
with the eVLA, (enhanced Very Large Array) but with substantially 
higher survey speed at frequencies below 1.4GHz. The Square 
Kilometre Array (SKA) is the next generation radio interferom- 
eter and will be a transformational instrument with an improve- 
ment in sensitivity compared to the eVLA of a factor of fifty with 
a 10^ increase in survey speed. At frequencies below 500MHz, 
aperture arrays will form the collector with dishes at frequencies 
above 1.4 GHz. In the mid frequency range (0.3 1.4 GHz) large 
survey speeds will be a chieved by using ape rture arrays or dishes 
with focal-plane arrays jSchilizzi et ahluOOTl) . Phase I of the tele- 
scope (with between 10 and 15 percent of the collecting area) is 
scheduled for operation in about 2015 with phased deployment of 
the rest of the telescope over the next few years and with comple- 
tion in 2021/22. Aperture arrays provide the most exciting technol- 
ogy for the mid-frequency range with the largest increase in survey 
speed provided they can mee t an appropriate performance to cost 
( [Alexander & F aulkner" 2009). The derived properties for th e SKA , 
based on lSchiiizzi et al.. (,2007i) and I Alexander & FaulkiieJ ( l2009b 
are summarised in Table[T] together with a comparison to the VLA. 

In this paper, we investigate the potential of these instruments. 



in particular the SKA, to measure cosmological evolution in cluster 
magnetic fields. We concentrate on Faraday rotation studies: the 
magnetic field along the line of sight rotates the polarisation angle 
by 

RM 



(1) 



where A is the observing wavelength and we included the cosmo- 
logical factor (z is the redshift of the rotating Faraday screen) to ac- 
count for the redshift of the radio emission traveling from the Fara- 
day screen to us. This factor ensures that for background sources 
with redshifts much higher than the targeted cluster, the contribu- 
tion to the rotation of the polarisation angle from the vicinity of 
the radio source will be negligible (unless that RM contribution is 
extremely high), since it occurs at much higher frequencies. RM is 
related to the thermal electron density, (cm^-'), and the line-of- 
sight magnetic field, B (flC), as: 



RM = 812 / «eS d/radm" 




(2) 



where the path-length dl is measured in kpc. 

The rotation measure has contributions from all along the 
line of sight, from the inter-galactic medium (IGM) in the vicin- 
ity of the source itself, from the intra-cluster medium (ICM) in 
the clusters along the line of sight, and from our own galaxy. The 
galactic contribution typically dominates at galactic latitudes below 
about \b\ < 20° . Above \b\ > 40°, t he galactic contribution becomes 
small er than |RM| < 30radm^^ jSimard-Normandin & Kronberd 
Il980l) . and can be as small as 5radm~^ jGuidetti et alj l2008h . 
Hence, high rotation measure from extragalactic sources at high 
galactic latitudes can most likely be attributed to Faraday rotation 
within the gas in galaxy clusters. 



1.1 Observations to date 

So far, the only individual cluster with a significant number of ob- 
served polarised members and background sources is the Coma 
cluster. Seven sources within 35 arc min display a clearly enhanced 
|RM | between 35 and 65 rad rcT^ jKim et al.|[l990l : iFeretti et all 
I1995I) . These sources, as well as control sources in the vicinity 
of the cluster, but avoiding sightlines through Coma's ICM, have 
been observed in many individual pointings with the VLA. For 
the cluster A514, six embedded radio sources have been observed 
in individual VLA pointings. These show a decline of iRMmaxI 
from 154 rad m^^ in the centre of the cluster to 15 rad m^^ on 
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Figure 1. Mass (top) and redshift (bottom) distribution of the adopted 
galaxy cluster model. See text for details of the model. 



the outskirts. Several other clusters, each with a f ew polarised 
sources, have been observed with similar results ('peretti et al.' 
ll999l : lTavIor~e t al. 2001; Govoni et al. 2006; Guidetti et al. 2008). 
Many clusters have only one or two sufficiently bright radio sources 
in or behind therrQ. Therefore stac king has been used to infer a gen- 
eral RM value for nearby clusters l lClarke et al.l200lh . A systematic 
increase in |RM| is found towards the cluster cores. 

First attempts have been made to infer the magnetic field struc- 
ture via the po wer-spectrum with F araday rotation studies. For the 
cluster A2255. TGovoni et al. measure the Faraday rotation 

of four embedded radio sources. For different assumptions about 
the magnetic field power spectrum they can extrapolate this data 
to predict the synchrotron emission of the radio halo. Comparison 
to observations of the radio halo then constrains the power law in- 
dex of the magnetic field's 3D-power spectrum to the range three 
to four, consistent with Kolmogorov turbulence, with a possible 
spatial variation. The model also predicts polarisation of the halo, 
which is indeed detected. 

This paper is structured as follows: First, we discuss the simu- 
lation method in Section|2] The Faraday rotation detection statistics 
and the implications for the measurement of the cosmological evo- 
lution of the magnetic field in clusters are presented in Section [3] 
We discuss our findings in Section|4l and summarize them in Sec- 
tion|5] Throughout, we use the following cosmological parameters; 
h = 0.7, Q.„r = 0.3, i^A = 0.7. 



See: |http : //www . mpa-garching ■ mpg . de/~kdolag/BCluster/ 1 



Polarisation model 




Figure 2. Distribution of the fractional polarisation for the radio source 
model. The dark blue triple dot-dashed Hne represents the fit from 
iTucci et alj )2003l) (not yet corrected for Faraday depolarisation). The solid 
black line represents our adaptation of this fit. This model is corrected 
for Faraday depolarisation by the polarisation dependent correction fac- 
tor shown as red dot-dashed curve (divided by 1000). The model finally 
adopted in this paper is shown as a green dashed line. The polarisation dis- 
tribution functions agree all at high fractional polarisations, were they are 
gauged against NVSS data. 

2 THE SIMULATIONS 

In order to assess the effect of cluster Faraday screens on the po- 
larisation of background radio sources, a relatively sophisticated 
model of clusters and source populations is required. 

2.1 Cluster model 

In order to model the cluster number within an observed patch 
of sky (here 10 deg^), we require a model for the cosmologically 
evolving cluster mass function. For this, we empl oy a cluster mass 
function whose multiplicity function (see e.g. ISheth & TormenI 
1 19991) is based on a moving barrier shape ansatz and depends on a 
minimal number of parameters to provide a suitable fit to numerical 
simulations. In particular, the shape of the used mass function has 
been derived from N-body simulations of large cosmological vol- 
umes and high resolutions following the method described in Jenk- 
ins et al. (2001). The above model gives the cluster mass function as 
a function of redshift. We adopt a lower mass cutoff of 3 X IO^^Mq 
and for a given realisation of the sky model derive cluster counts as 
a function of redshift, making allowance for cosmic variance. Our 
resulting galaxy cluster catalogue contains 775 objects, with red- 
shifts up to --^ 2.3. The properties of our generated galaxy cluster 
sample are illustrated in Fig.[T] 

In our model, clusters are distributed randomly and radio 
sources are either background or foreground objects with respect 
to a given cluster. Every object is placed on a 3D grid, covering 
\/T0 X \/T0 square degree and the redshift space from zero to eight. 
The spatial resolution is 6000 x 6000 cells, and redshift space is 
covered by 240 cells. This gives an effective resolution of 1.9 arc- 
sec and dz = 0.033 in the third dimension. 

2.2 Rotation measure model 

The observational data on the st ructure of magnetic fields in clus- 
ters is currently sparse. Based on IClarke et al.H200lh . we therefore 
employ a simple model for the rotation measure, with a peak value 
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Rotation measure map 
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Figure 3. Top: Derived distribution of cosmological rotation measures on tlie sky due to galaxy clusters. Bottom: Recovered RM-sky due to background radio 
sources for a 100 h SKA-AA observation. To enhance contrast, the resolution has been reduced by a factor often compared to the figure on the top. The real 
field of view as planned to date would be 25 times larger. 
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(justified in Section [Z9l below) drawn from a uniform random dis- 
tribution between RMq = ±5000 rad m^^, and a decline according 
to a j3 -model: 

RM(r) =RMo(l + rV'-o)"^*"''^- (3) 

We use /jRM = 1 (see below) and ro = 0.0636 r2oo 

teiviano & Saluccil l2006h . where r2()o is the radius within 
which there is an over-density of a factor 200 compared to the 
mean density of the universe at the redshift of the cluster. r2oo is 
roughly equal to the virial radius of the cluster. This results in core 
radii between 20 and 110 kpc, which corresponds to angular sizes 
bet ween a few arcseconds to h alf an arcminute. 

iBiviano & Saluccil ( |2006|) give a value of /3 = 2/3 for the gas 
density distribution in clusters of galaxies. The temperature distri- 
bution in clusters of galaxies is typically not far from isothermal 
(except in the centres of cool core clusters). Therefore, the pressure 
follows a very similar radial profile as that of the density. Equipar- 
tition between the magnetic and the thermal pressure is probably a 
reasonable assumption for the magnetic field distribution in clus- 
ters, and consistent with some RM-data (e.g. iGovoni et al]|2006l : 
iGuidetti et al ]|2008). Consequently, we adopt a value of jS^ = 1 /3 
for the radial distribution of the magnetic field strength along the 
line-of-sight. The line-of-sight integral then results in the /3-model 
for the rotation measure described in equation l[3} above, with 
|8rm = 1 • We calculate the rotation measure contribution out to the 
virial radius. 

We show in Fig. |3] (top) the resulting distribution of rotation 
measures due to the simulated galaxy clusters. Much of the sky is 
covered with galaxy clusters, however, many of the systems have 
small rotation measures. About 12 per cent of the sky is covered 
with regions of at least 10 rad m^^. The fraction reduces by a factor 
of about two if one considers rotation measures above 20 rad m^^. 
In our chosen cosmology, only 34 per cent of all clusters are lo- 
cated beyond a redshift of one. These are obviously smaller on av- 
erage than the lower redshift population. Hence, essentially all of 
the area coverage happens up to a redshift of unity (Fig.|4ll. There- 
fore, up to this redshift, structure determination should be much 
easier than beyond. 

2.3 Radio source model 

To obtain realistic models of the radio source populations, we use 
the publically available SKADS (Square Kilometre Array Euro- 
pean Design S tudy) synthetic radio s ource catalogue {SKADS simu- 
lated skies, 5^: IWilman et alj200 j^l). The catalogue is derived from 
semi-empirical simulations of different populations of radio contin- 
uum sources and extends down to faint flux limits in order to allow 
observation simulations of high-sensitive future radio telescopes, 
such as the SKA. The source populations included in the catalogue 
are radio-loud (FR I and FR II) and radio-quiet AGNs as well as 
radio emitting starburst and normal disk galaxies. The luminosity 
function corresponding to each class of sources is derived empiri- 
cally from available observations. The fits are extrapolated to low 
luminosities and assume redshift evolutions to assure the simula- 
tions are complete down to the instrumentally expected (faint) flux 
limits. Our version of the catalogue covers ten square degrees of 
sky out to a r edshift of eigh t with a 1.4 GHz flux limit of 1 /xJy in 
Stokes/. See lWilman et al.l ( l2008h for a detailed description of the 
radio source simulations. 

^ See: |http : // s-cubed . phy s Ics . ox . ac . uk7] 



2.4 Polarisation model 

The degree of polarisation has to be determined for individual 
sources. We base our prediction of the polarisation of these radio 
sources on the statistics of the NRAO VLA Sky Survey (NVSS) 
catalogue (Condon et al. 1998). Tucci et al. (2003) have analysed 
polarisation statistics from the NVSS and other sources in detail. 
Depending on spectral shape and source flux, they derive median 
fractional polarisations between 1.12 and 1.77 per cent for NVSS 
sources. For the NVSS (1.4 GHz), the probability distribution of 
the fractional polarisation declines monotonically, but for higher 
observing frequencies, a turnover may be seen towards lower frac- 
tional polarisation. This is due to a systemat ic increase of the f rac- 
tional polarisation with observing frequencv. lTucci et al.l ( l2003l) ex- 
plain these findings by Faraday depolarisation. Therefore, we ex- 
pect that the low frequency distribution function also turns over at 
a not yet measured fractional polarisation. We therefore fit a log- 
normal function to the measured part of the fractional polarisation 
distribution function at 1.4 GHz. This results in a med ian value of 
1.64 p er cent and a standard deviation of 0.7 per cent. iTucci et al] 
( I2OO3I) believe that some of the NVSS sources are depolarised by 
differential Faraday rotation in intervening galaxy clusters. We ar- 
gue that this would be much less so for SKA observations: The 
main reason is that with lower flux limit the radio source popula- 
tion changes. Extended jet sources dominate at high flux densities, 
and star forming galaxies at increasing distance and decreasing an- 
gular size at lower flux densities. For high redshift galaxy clusters, 
the depolarisation frequency is reduced by a factor (1 Ex- 
tended sources behind nearby clusters will be much better resolved 
than by the NVSS (45 arcsec). For a turbulent power spectrum, the 
RM variations should have less power on smaller scales, h ence the 
depolarisation frequency will be lower jKrause et alj2007l) . 

For the se reasons, we adopt the de-depolarisation model from 
iTucci et all (2003) to derive the intrinsic polarisation H from the 
potentially Faraday depolarised polarisation HpD, that is drawn 
from a random distribution as described above: 

n = nFD/(nFD,VGHz)//(nFD,l-4) (4) 

f{p.y) = 721og(0.5v'''* + 0.75) exp(-3.2/-^^ +0.8). 

Here, Vghz 's a reference frequency in GHz. We use 10 GHz as 
the reference frequency, assuming that Faraday depolarisation for 
NVSS sources mainly occurs below 10 GHz. The correction in- 
creases the fractional polarisation by a factor of about two, for those 
sources with low polarisation at 1 . 4 GHz, as shown in Fig. |2] This 
effect has been measured already: iTavlor et alJ (l2007l) measure an 
increase of the fractional polarisation downwards of 100 mjy of a 
factor of three compared to the high flux density sources, in agree- 
ment with the model we propose. 



2.5 Redshift distribution 

For decreasing flux limit, the radio source counts are dominated 
by sources of increasing redshift. We show the redshift distribution 
for different (polarised) flux limited samples in Fig. [5] For a po- 
larised flux-density of 100 ^Jy, roughly the level of current studies 
dTavlor et al.ll2007r) . the median redshift is 0.6. It rises to 1.3 for a 
polarised flux limit of 0.1 ^Jy, appropriate for a 100 h SKA obser- 
vation. Even for the highest redshifts, we still predict thousands of 
background sources per square degree for this flux limit. 
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2.6 Required signal-to-noise for RM measurements 



2.7 Rotation measure grid 



A distinctive feature of the next generation of radio telescopes is 
that a large bandwidth and high spectral resolution will be avail- 
able simultaneously. Rotation measure synthesis is an established 
technique to analyse RM data. However, for the analysis of the RM- 
grid, we wish to extract information from the faintest detectable 
sources. For the faintest sources we will extract a single value of the 
rotation measure. Here, we consider the minimum signal-to-noise 
on a point source required for this. Faraday rotation leads to sinu- 
soidal variation with respect to across the band for the Stokes 
parameters Q and U, modulated by a function giving the fractional 
polarisation as a function of wavelength and the intrinsic spectrum 
of the source. We consider normalised Stokes parameters: 



u' 



-1/2 



(5) 



cos(2PA + 2RMA2) 
sin(2PA + 2RMA2) 



where Q and U are the ordinary Stokes parameters. To assess the re- 
quired signal-to-noise ratio (SNR), required of a background source 
to measure a given foreground RM, we perform Monte-Carlo sim- 
ulations. Our simulated data were for a background source with 
a spectral index of 0.7 - we simulate the observing band, taking 
a constant channel width of Av = 1 MHz. We perform a com- 
bined X^'^^ to the normalised Stokes parameters, and use the flux 
at 1 .4 GHz as a reference point for the SNR. For infinite SNR, the 
minimisation of would be straight forward. For decreasing SNR, 
additional minima appear in x^ as a function of RM. For a SNR of 
10 (100) the distance between these minima is of order 1 rad m^^ 
(100 rad m^^). A dense grid of starting guesses is therefore essen- 
tial, for a gradient search method to find the global minimum. The 
actual values of X^ differ only by a few percent, with the excep- 
tion of the x^ for the true RM, and its negative. We find that x^ for 
the true RM and its negative, are significantly lower than all other 
RMs (five times the standard deviation for 400 starting guesses) for 
a signal to noise of four and higher 

The sign ambiguity is due to the symmetry properties of 
the trigonometric functions. If a set (PA,RM), where PA is the 
high frequency polarisation angle (compare equation (O), are the 
correct parameters to describe both functions, 2'(PA,RM) and 
(/'(PA.RM), then (;r + PA, -RM) will be another solution regard- 
ing Q' and (— PA,— RM) respectively for U': 



g'(PA,RM) 
i7'(PA,RM) 



2'(;r + PA,-RM) 
(/'(-PA,-RM) 



To fix the sign at low SNR, for each set of fitted parameters for 
Q' we also performed two fits for U' with starting guesses of 
(;:: — PA,RM) and (PA,RM), respectively, allowing only RM to 
vary, i.e. keeping the high frequency polarisation angle fixed. If 
the fit with (PA,RM) leads to the lowest minimum of x^^ then we 
have found the correct RM. However, if the fit with (;r — PA.RM) 
produces the lowest minimum of x^, then — RM will be the cor- 
rect rotation measure. The procedure is then repeated with U' and 
2' exchanged. We may thus produce a combined x^ for RM and 
— RM. Selecting on the combined x^^ we are able to recover the 
correct sign. 

For RMs between 3 and 3000 rad m^^, we find a minimum 
SNR between four and ten. This is true for the observing bands 
from 800 MHz to 1050 MHz as well as 300 MHz to 1000 MHz. In 
the following, we assume that a RM measurement may be done for 
any source detected at an SNR of ten in polarised flux. 



For observations of nearby clusters, as well as the galaxy, it is useful 
to compute the root-area-average distance between polarised back- 
ground sources. Fig.[6]shows the average grid spacing over observ- 
ing time for the adopted SKA and SKA phase 1 sensitivities. We 
use the point source sensitivity limit according to lWrobel & Walked 
( Il999l) : 



kBTsy 



\/2dv tAe 



(6) 



where, is Boltzmann's constant, Tsys is the telescope's system 
temperature, A the total effective area, dv the bandwidth, T the 
observing time and e the overall efficiency. Adopting a SNR of 
ten, we count the sources above the sensitivity limit for a particular 
telescope and observing time. The total source counts are roughly 
inversely proportional to the polarised flux, which is an effect of 
the adding up of the contributions from different source types. So, 
while the galaxies constitute the majority of the sources at low 
fluxes, other source types still contribute enough to flatten the over- 
all distribution considerably. 



2.8 Expected RM values 

For each background source, we calculate the observed rotation 
measure by integrating through the simulation volume from the ob- 
served source to the observer We show the cumulative histogram of 
all the rotation measures produced in this way in Fig.|7] We find that 
the majority of rotation measures are in the range 10-100 rad m^^, 
due to the fact that most of the area comes from the outskirts of the 
clusters. The highest rotation measures come of course from the 
cluster centres. For a given frequency band and rotation measure 
range, a particular SNR is required. As shown above in Section|2j6] 
both of the frequency ranges, 300 MHz to 1000 MHz (aperture ar- 
ray) and 800 MHz to 1050 MHz (dish array), should be suitable for 
the great majority of expected rotation measures. 



2.9 Core rotation measures 

For each radio source within three virial radii of a given cluster cen- 
tre, we also calculate the respective impact parameter, i.e. the pro- 
jected distance to the cluster centre. We show these rotation mea- 
sures against the impact parameter, for low redshift (z < 0.5), high 
mass (M > lO'^M^) clusters and a pol arised flux limit of I /iJy 
in Fig. [S] These cuts were employed in IClarke et alJ tOOA and 
their data is also shown in Fig. [S] But our model and the obser- 
vational data show only RM values up to |RM| < 300 rad m^^, 
whereas our model has peak values of 5000 rad m^^ in the clus- 
ter centres. The reason is that at the high flux cut adopted here, 
the RM-grid is still very sparse, and hardly any source is located 
within the core radius. If we employ a lower central RM, with the 
same cuts, the simulated RM spread on the scale of the virial radius 
would be too small compared to the observed data. As can be seen, 
the central rotation measures of about 5000 rad m^^ that are used 
in the simulated data here are necessary to reproduce the observed 
spread of rotation measures. For distances greater than 1000 kpc 
and random sightlines, we enco unter high RMs from neighbouring 
clusters. These are not present in lClarke et alj ( [200lh . as they select 
their control sample avoiding known X-ray clusters. 
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Figure 4. Fraction of sky covered with RM> 10 rad m^^ (upper curve) and 
RM> 20 rad m"- (lower curve). 
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Figure 5. Distribution of radio sources against redshift for different po- 
larised flux limits. From top to bottom, the curves represent: 0.1, 1, 10, and 
100 fiJy. The median redshift is 0.6, 0.6, 0.8, and 1.3. 
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Figure 6. RM-grid spacing for different telescopes and a SNR of ten. 
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Figure 7. Cumulative histogram of all the rotation measures we expect to 
find within one square degree for a polarised flux limit of 1 fl]y (lower red 
curve) and 0.1 jlJy (upper black curve). In any case, we expect most of 
the measurements to yield rotation measures of order 10 to one hundred 
rad m^^. 



3 RESULTS 

3.1 Detection statistics 

We show the detection statistics for the SKA aperture array and the 
phase 1 SKA in Fig. [10] and Table |2l The histograms demonstrate 
the progress that can be expected from the SKA, especially the full 
aperture array. 

For a one hour phase 1 dish blind survey, we predict to find 
30, 214, and 148 clusters in the redshift ranges 0-0.2, 0.2-0.5, 
0.5-1, respectively, with more than one background RM measure- 
ment, each. We also predict a few detections at higher redshift. This 
would allow stacking experiments for different redshift bins out to 
a redshift of at least one. Also, one could determine the decline of 
the average RM with cluster radius better than presently possible. 
This could answer the question weather the magnetic field energy 
density follows the thermal energy density. We also would have 
some chance to find about a couple of nearby clusters with 30-100 
RMs, allowing an improved determination of the field structure. 30 
back ground sources should be enough to cons train the power spec- 
trum jGovoni et alj|200d iGuidetti et al1l2008h . Since the positions 
of the big nearby clusters are known, a targeted survey would cer- 
tainly yield about this number of RMs. 

A 100 hour phase 1 observation would result in 26 clusters 
at z < 0.5 with more than 30 background sources, and a couple 
of these having more than 300 background sources. This would 
allow to study the field structure in some detail. Up to a redshift 
of unity, nearly every cluster in the field of view should have at 
least one background source. Even beyond z = 1, we predict to find 
490 clusters with one to several background sources. This would 
allow the extension of the stacking method to this redshift. 

The full SKA aperture array will finally allow to measure the 
structure of the magnetic field in great detail: Even for a 1 hour 
observation, each low redshift cluster will at least have 10 back- 
ground sources 25 of which over three hundred. We predict to find 
325 clusters in total with more than 30 background sources, dis- 
tributed over a redshift range up to 0.5, promising good statistics 
of the variance of the power spectrum from cluster to cluster. More 
than 10,000 clusters would be detected with at least one RM mea- 
surement beyond a redshift of 0.5. 

For a 100 hour pointing, we will be able to constrain the power 
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impact parameter b / kpc 

Figure 8. Rotation measure against impact parameter. The impact parame- 
ter is the projected distance of a radio source to a given cluster centre. Low 
redshift {z < 0.5), high mass (M > 10''* Mr ) clusters have bee n selected 
in order to compare to the o bserved data of [ Clarke et al] jlQOlh . The po- 
larised flux l i mit is 1 jlJy, as in IClarkeet alj)200lh . Observational data from 
IClarke et al.l j200lh (cluster sample only) is provided as red crosses. For 
distances greater than 1000 kpc, the sightlines to the radio sources pierce 
neighbouring clusters. To get the observed RM spread at the observed im- 
pact parameters, central RMs of about 5000 rad m"' are necessary. 
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Figure 9. Maximum rotation measure seen against any given cluster in our 
sample against cluster redshift for a polarised flux limit of 0. 1 /iJy. The red 
lines denote the median for sub-samples with corresponding redshifts (0- 
0.5,0.5-1, 1-1.5, 1.5-2, > 2). The purple lines denote the measured rotation 
measures if the true rotation measures would decline with redshift as ( 1 + 
z)^', the orange line displays the medians for an increase proportional to 
(l+z). 

spectrum for many clusters in ttie field of view up to a redshift 
of one (> 5,000 clusters with more than 30 background sources). 
Up to a redshift of one, we expect > 1 , 000 clusters with at least 
100 RMs. Only for such clusters, we may expect to resolve the 
core in some cases, i.e. to have at least one background source 
within the cluster's core radius. A 100 hour pointing would corre- 
spond roughly to the sensitivity limit used in Fig.|9l where we show 
the decrease of the maximum rotation measure in any given cluster 
with redshift due to the decreased likelihood to get a background 
source within the core radius. 25 clusters should have more than 
1000 background sources, and hence allow a very detailed deter- 
mination of the magnetic field structure, even in the cluster cores. 
Due to the limited field of view, the full SKA dish array would 



not nearly reach the detection rates of the aperture array. Even a 
100 hour blind survey would only yield 24 clusters with more than 
30 background sources. 

3.2 Cosmological evolution 

Clusters at higher redshift cover less area of sky (up to a redshift 
of 1.6). Also, there are fewer sources behind a unit area at higher 
redshift. Therefore, the probability that a cluster centre has a back- 
ground object decreases as the redshift increases. We plot the max- 
imum rotation measure that we detect for a polarised flux limit 
of 0.1 fl]y against any given cluster versus the cluster redshift in 
Fig.|9] Out to a redshift of unity, there is still a good chance of find- 
ing sightlines through the cluster core. However, the median maxi- 
mum rotation measure decreases by almost a factor of four for clus- 
ters with redshift 0.5-1 compared to ones at 0-0.5. The decrease in 
the median of the measured maximum rotation measure is almost 
exponential with redshift, with an 1 / e-decay length in z close to 
0.4. We can check the detectability of possible trends with redshift 
by multiplying the rotation measures with (1 where we have 
chosen n to be ±1 as example. The change in the median rotation 
measure is also indicated in Fig. |9] The decay length changes to 
0.3 and 0.5, respectively. If we compare the redshift intervals 0-0.5 
and 0.5-1, we predict a decrease of the median rotation measure 
by factors of about 3,4 and 5, respectively (compare with Table[3] 
which is normalised to the low redshift value). 

The statistical error on the median maximum rotation measure 
is given by: 

Om = l.253a/VN, 

where <y is the standard deviation and A' is the number of sources in 
the respective bin. We have calculated the median maximum RM 
for any given cluster and the statistical errors for a 100 h SKA AA 
observation for different values of n. This is shown in Table[3] Dif- 
ferences in the exponent n of 0.3 could be distinguished. 



4 DISCUSSION 

We assess the capability of the SKA and its precursor, phase 1 SKA, 
to detect rotation measures and determine the evolution of magnetic 
fields in galaxy clusters. Our cluster and radio source simulations 
are based on standard cosmological methods and extrapolation of 
available observational data. We use a fractional polarisation model 
based on NVSS data, moderately corrected for Faraday depolarisa- 
tion, at low fractional polarisations, only. 

We find that about 10 per cent of the whole sky should be 
covered with rotation measures greater than 10 rad m^^. If the 
cosmology is known, this may serve as an additional constraint to 
check for the proper removal of the Milky Way's RM foreground. 
Large clusters can have a RM contribution from background clus- 
ters. However, the contamination occurs mainly on the outskirts of 
these clusters. 

Our simulations show that core RMs of about 5000 rad m^^ 
are necessary to explain the currently known stacking data, which 
shows about 200 rad m^^ on the scale of the virial radius. RMs 
against central cluster radio sources often reach values comparable 
to our core RMs ( Carilli & Taylor 2002). This suggests that our 
assumption for the core RMs is realistic, but also that the functional 
form of the radial decline we assume is not too far from reality. 

How would the derived statistics change if one were to change 
the radial RM-profile? The adopted profile has high core values 



Cluster magnetic fields with the SKA 9 



Numbers of clusters per pointing with differing numbers of bacl^ground sources 
SKA Phase 1, Ihr 



10000 




I 0<z<0.2 






I 0.2<z<0.5 






I 0.5<z<1.0 






I z>1.0 



^ 100 - 



O 10 - 



Number of background sources per cluster 



Numbers of clusters per pointing with differing numbers of background sources 
SKA Phase 1, 100hr 



I 


0<z<0.2 


I 


0.2<z<0.5 


■ 


0.5<z< 1.0 


I 


z> 1.0 



Number of background sources per cluster 



Numbers of clusters per pointing with differing numbers of background sources 
Full SKA Dishes, Ihr 



lli 



Numbers of clusters per pointing with differing numbers of background sources 
Full SKA Dishes lOOhr 



I < z < 0.2 




10000 




I 0<z<0.2 


I 0.2<z<0.5 








I 0.2<z<0.5 


I 0.5<z<1.0 








I 0.5<z<1.0 


I z > 1 .0 








I z>1.0 



Number of background sources per cluster 



Number of background sources per cluster 



Numbers of clusters per pointing with differing numbers of background sources 
Full SKA AA, Ihr 



^ 100- 

Q. 

B 
"55 

d 10 - 



I 


< z < 0.2 


I 


0.2<z<0.5 


I 


0.5<z< 1.0 


I 


z>1.0 



Number of background sources per cluster 



Numbers of clusters per pointing witti differing numbers of bacl<ground sources 
Fuii SKA AA, lOOhr 



O 10 



I 


0<z<0.2 


I 


0.2<z<0.5 


I 


0.5<z< 1.0 


I 


z> 1.0 



Number of background sources per cluster 



Figure 10. Detection statistics for the SKA aperture array (bottom row), the full SKA dish array (middle row) and the phase 1 SKA dish array (top row), 
telescope specs from Table[T] for an observing time of 1 hour (left) and 100 hours (right), respectively. We plot the number of background sources per cluster 
on the horizontal axis, and the cluster counts for the particular bins on the vertical axis. The huge number counts found with the full SKA aperture an'ay is a 
merit of its huge field of view (250 square degrees) combined with its high sensitivity. The same information is provided in tabular form in Table|2] 
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Table 2. Cluster counts with a given number of background sources. 



Redshift 






Number of polarised background sources per cluster 






interval 
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0-0.2 








8 


20 





2 











0.2-0.5 


54 


160 


54 




















0.5-1 


568 


148 























> 1 


522 


10 























Phase 1 Dishes, 100 hour pointing 


0-O.2 














8 


20 





2 





0.2-0.5 








4 


162 


102 


4 











0.5-1 


2 


86 


444 


182 


2 














> 1 


42 


370 


118 


2 

















Full SKA dishes, 1 hour pointing 


0-0.2 














1.6 


1.2 





0.2 





0.2-0.5 





0.2 


8.6 


14.6 


3.8 














0.5-1 


1 


32 


35.8 


2.8 

















> 1 


15.2 


35.4 


2.6 




















Full SKA dishes, 100 hour pointing 


0-0.2 




















2.4 


0.4 


0.2 


0.2-0.5 














8.6 


16.8 


1.8 








0.5-1 











20.2 


49.2 


2.2 











> 1 





0.2 


10.2 


39.4 


3.4 














Full SKA AA, 1 hour pointing 


0-0.2 














100 


250 





25 





0.2-0.5 








50 


2150 


1150 


50 











0.5-1 


25 


1175 


5550 


2175 


25 














> 1 


575 


4750 


1300 


25 

















Full SKA AA, 100 hour pointing 


0-0.2 




















150 


200 


25 


0.2-0.5 














125 


2475 


800 








0.5-1 











350 


6925 


1675 











> 1 








150 


4575 


1925 















with a steep decline. With this distribution we have shown to match 
the observations of nearby clusters. However, one might argue that 
these nearby clusters may not be representative of whole cluster 
population. One might imagine that a less centrally peaked flatter 
distribution, like for example in the Coma cluster, might be more 
typical. We have shown in section lZ6l above. that the optimum RM- 
range is between a few and a several 100 rad roT^. Such a scenario 
would therefore lead to an increase of the fraction of sky in our 
best range. Therefore, the RM-distribution we have chosen is, if 
anything, rather pessimistic. 

We show that with a phase 1 SKA, as defined in Table[T] good 
stacking experiments should be possible, determining the average 
magnetic field profile locally and also out to a redshift of about one. 
The field structure should be measurable for some nearby clusters. 
If there would be a choice necessary between a deep survey and 
100 one hour pointings, one should probably do the deep survey, 
since it already allows some field structure determination in nearby 
clusters, and promises to measure RMs for clusters with z > 1. 

With the full SKA AA aperture array, both a shallow 100 one 
hour pointing survey, and a deep survey promise very interesting 
results. A shallow 100 times one hour survey should detect RMs 



from about 10^ clusters, with complete redshift coverage. A deep 
100 h one field survey will allow structure determination in un- 
precedented detail out to a redshift of one: over 1000 clusters will 
have more than 100 RM measurements, 25 of which over 1000 
RMs. The result of a 100 h deep survey is also shown graphically 
in Fig.|3](bottom). To some extent, the two approaches are comple- 
mentary: the 100 h deep survey would offer a complete RM survey 
for all structures above 3 x IO'^Mq, and thereby exclude possible 
selection biases which could arise if only a fraction of the clusters 
is detected as in the shallow survey. The shallow survey would be 
able to detect large scale bias, and also would provide better overall 
statistics on the cosmological field evolution, if details within the 
clusters are neglected. 

Without phased array feeds on the majority of the dishes, as 
assumed here, the detection rates remain rather small for the full 
SKA dish array. However, having similar sensitivity than the aper- 
ture array, but coverage at higher frequencies, it would be a valuable 
complement for high RM regions (cluster centres), where the fre- 
quencies accessible to the aperture array would be depolarised.The 
RM grid density would be similar to the aperture array. 

With the deep survey, we could measure the cosmic evolution 
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Table 3. Median maximum rotation measure in any given cluster against redshift for different 
assumptions about the intrinsic evolution, where n is the exponent for a systematic change 
with redshift. The statistical errors are applicable for a 100 hour SKA AA observation. The 
values in each row have been normalised to the lowest redshift bin. 



n 


z = 0-0.5 


z = 0.5-1 


z= 1-1.5 


z= 1.5-2 


z>2 


1 


1± 0.015 


0.346 ±0.010 


0.168±0.012 


0.066 ±0.009 


0.024 ±0.058 


0.5 


1 ±0.015 


0.298 ± 0.009 


0.128 ±0.009 


0.047 ± 0.007 


0.016±0.038 


0.2 


1 ±0.015 


0.276 ±0.008 


0.109 ±0.007 


0.039 ±0.005 


0.012 ±0.029 





1 ±0.015 


0.259 ±0.008 


0.099 ±0.007 


0.034 ±0.005 


0.010 ±0.024 



of the RMs in clusters. If RM is proportional to (1 + z)", we could 
measure n to an accuracy of An = 0.3. In order to determine the 
evolution of the magnetic field from this, one needs additional in- 
formation on the electron densities. If we would adopt an accuracy 
for n of An = 0.3, and assume a similar behaviour and accuracy for 
the ele ctron densities, to be measured via the Sunyaev-Zel'dovich 
effect dSehgal et alJlOOSi) . we would end up with an accuracy for a 
power law index for the magnetic field (6 <>= (1 +z)"') of 0.4. This 
should be regarded as a conservative estimate, since it only takes 
into account information from the deep survey. Additional infor- 
mation from a shallow survey should improve the accuracy for the 
low redshift bins. Also, we have only used the information on the 
maximum RM per cluster. Taking into account the full data, will 
also improve the statistics. 

For a polarised flux limit of 0.1 /xJy, as appropriate for a 100 
hour SKA pointing and a SNR of ten, we predict a few 10,000 
radio sources per square degree. These sources have a broad red- 
shift distribution, with a median redshift of about one, rising with 
decreasing flux limit. We can hence expect not only to find the re- 
ported numbers of background sources, but also a sizable number 
of foreground sources, of order a few percent of the background 
sources, rising of course with cluster redshift. An accurate deter- 
mination of the foreground RM opens up the possibility to detect 
a possible small mean cluster RM which would point to a super- 
cluster scale field connecting the cosmic web. This would be the 
first detection of such a field. 



5 CONCLUSIONS 

We have modelled the distribution of rotation measures for the 
SKA. We calculate the average distance between rotation mea- 
sures to be between 6 arcminutes for a one hour observation with 
the early SKA and 36 arcseconds for a hundred hour pointing 
with the full SKA. We expect to find rotation measures mostly 
up to 100 rad roT^, with the cluster centres reaching up to several 
thousand. The planned SKA mid frequency aperture array (300- 
1000 MHz) as well as the lowest band of the dish array would be 
well suited for the great majority of sources. High RM cluster cen- 
tres would require targeted high frequency follow up observations. 
We find that current Faraday rotation studies only represent the out- 
skirts of galaxy clusters. The average RM should increase steadily 
towards the core to reach typically several thousand rad nT^ . 

A phase 1 SKA would already improve the statistics for cur- 
rent stacking experiments considerably, and make this experiment 
viable out to redshifts z > 1. The full SKA aperture array would 
detect over a million clusters with at least one background source 
each in a shallow 100 hour survey, and allow detailed field structure 
determination (> 1000 clusters with more than 100 background 
sources each) with a deep survey. If the cosmological evolution of 



the rotation measures is proportional to (1 + the SKA would be 
able to measure n to an accuracy of 0.3. Compared to the few RMs 
known for a few nearby clusters today, this will revolutionise our 
knowledge of rotation measures in galaxy clusters. Provided the 
electron densities can be measured via the Sunyaev-Zerdovich ef- 
fect at high redshift, we can expect to follow the build-up of cosmic 
magnetism with the SKA. 
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